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Orbital symmetry of a triplet pairing in a heavy Fermion 
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The orbital symmetry of the superconducting order parameter in UPts is identified by evalu- 
ating the directionally dependent thermalconductivity and ultrasound attenuation in the clean 
limit and compared with the existing data for both basal plane and the c-axis of a hexagonal 
crystal. The resulting two component orbital part expressed by (Aa;(k), Ay(k)) is combined with 
the previously determined triplet spin part, leading to the order parameter of either the non- 
unitary bipolar state of the form: d(k) = bAa;(k) -I- ijAj,(k) or the unitary planar state of the 
form: d(k) = bAi:(k) -f jAy(k) where b ± j = c, or a with the hexagonal unit vectors a, b and c. 
The d vector is rotatable in the plane spanned by a and c perpendicular to b under weak applied 
c-axis field because of the weak spin orbit coupling. Experiments are proposed to distinguish 
between the equally possible these states. 

KEYWORDS: heavy Fermion superconductor, UPta, triplet pairing, orbital symmetry, thermalconductivity, ultra- 
sound attenuation, bipolar state, planar state 
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§1. Introduction 

Much attention has been focused on various uncon- 
ventional superconductors in the recent years, includ- 
ing high Tc cuprates, ruthenate, low dimensional organic 
conductors and heavy Fermion materials. Extensive ex- 
perimental and theoretical studies on UPta, a typical 
heavy Fermion superconductor, have been devoted to in 
order to elucidate tte pairing symmetry realised in this 
particular materiallil^ UPts discovered in 1983B is rather 
special in this class of materials because it is found in 
lOSaS-* that the superconducting transition temperature 
is split (Tci ~0.58K and Tc2 -0.53K), self-cvidcntly sug- 
gesting a rich internal structure of the Cooper pair, def- 
initely not a simple conventional s-wave pairing. There 
are at least three phases identified: A (high T and low 
H), B (low T and low H) and C (high T and high H) 
phases. In spite of these facts and the following much 
efforts, however, it was a difficult task to identify the 
exact pairing symmetry. ^ 

Vast majority of experimental worksEl done on UPts 
are concerned with the energy gap structure in the su- 
perconducting state by measuring the T-dependences of 
thermodynamic and transport quantities whose power 
law behaviors at low T's point to the existence of at least 
a line node at the equator of the Fermi surface. Precise 
gap structure in addition to that line node, however, re- 
mains still to be worked out. At this stage oLthe reSj 
there weifi_tkree major scenarios by Joynt^* Sauli 
and ours&BllJ' who try to consistently explain the above 
Tc-splitting phcnomenon,|-Qther than the so-called acci- 
dental degeneracy model.tll' 
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In general, the pairing symmetry consists of the or- 
bital part and the spin part in the pairing function: 
AQ_^(k) — i{a ■ d(k)(T2)Q,/3- The former determines the 
gap topology which the above measurements probe. The 
latter characterizes the spin structure of the Cooper pair, 
which is related to the spin susceptibility of the system, is 
usually probed by the Knight shift experiment of NMR. 

A long waited Knight shift experiment by using high 
qualiiM jSingle crystals has been finally performed by Tou 
et al,ll3'll3) revealing the rich details of the spin structure 
of the pairing function in UPts: (1) The spin suscepti- 
bility Xi for the i-direction decreases below T < Tc2 for 
i — b and c-axes of the hexagonal crystal at low fields, 
corresponding to the B phase. (2) The parity of the pair- 
ing function is odd and the spin triplet state is realized. 
(3) The low temperature and low field state B phase 
is characterized by having the two-component d- vector. 
By increasing H (|| c) one of the d- vector component 
rotates to lower the Zeeman energy around 2kG, imply- 
ing that the effective spin-orbit coupling which is to pin 
the d- vector to JJia^ crystal lattice is weak. Thes£L.ex- 
perimental factsllSllj-' fit to our predicted-scenario p □EH' 
negating the £2^ scenario due to SaulsoQ' who assumes 
the onc-componcnt d-vector fixed to the c-axis by strong 
spin-orbit coupling and also ruling out the Eig scenario 
due to Joyntu' who predicts the spin susceptibility to 
decrease for all iJ-orientatioiia- This is not observed by 
Tou et al's NMR experimentlij' 

So far we have not so seriously considered the orbital 
structure of the pairing function partly because without 
identifying its spin structure it was difficult to narrow 
down the precise orbital symmetry among vast possible 
candidates. Now since its spin structure, we believe, has 
fairly well determined, we must select out the orbital 
part to complete our task. 
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For this purpose we calculate the T-dependences of 
the thermalconductivity and ultrasound attenuation co- 
efficient, each for major crystalline orientations of the 
c-axis and the basal plane for the two candidate orbital 
states whose spin structures are identical, namely, the 
non-unitary bipolar state and the unitary planar state. 
By adjusting several material parameters we try to fit 
the data of the above quantities in a consistent manner. 

Arrangement of the paper is as follows: We first in- 
troduce the two candidate orbital states, both of which 
can equally explain the above Knight shift experiment 
and thus logically same in this respect. In next Section 
we calculate the thermalconductivity and ultrasound at- 
tenuation coefficient by using these states. In the Sec. 4 
we explain the phase transitions in terms of these states. 
The final Section is devoted to discussions and conclu- 
sions which summarize our present status for the identi- 
ficatim|i_nL|the pairing symmetry done over the past ten 
years.l'Blia) 

§2. Non-unitary bipolar state and unitary pla- 
nar state 

2. 1 Choice of the order parameters 

The inert p-wave states with odd-parity under the 
hexagonal symmetry are enumerated group-theoretically 
when the spin-orbit coupling is weak,t3' which is the 
present case. An extension to more higher order an- 
gular momentum state is trivial, simply replacing the 
irreducible two-dimensional representation whose basis 
function is given by (fca,, ky) for the p-wave case in the 
cubic symmetry to the corresponding two-dimensional 
Eiu and 'Ei2u representations in the hexagonal symme- 
try characterized by {\xi}^)^ Aj^(k)). The enumerated six 
states are listed in Table 1. 



Table I. The enumerated states belonging to the two-dimensional 
representations characterized by {\x(}i.), Ay(k)) where r = iaa2- 
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Among these states listed in Table 1 several require- 
ments in light of the UPta problem should be met to se- 
lect. The decisive one is that the spin component must 
be two according to the above Knight shift experiments. 
PreviousbL we advocated a /3 type state or a 7 type 
stateJa'Q'H but these turn out not to be adequate by the 
foUoadng two reasons: (A) The /LtSR experiment by Luke 
et alliS' claims that the spontaneous moment appears be- 
low Tc2 , suggesting the time reversal symmetry breaking 
state realized as the B phase. The non-unitary /? state 
is consistent with this, but not the planar and bipolar. 



However, no later experiments reproduce their result. tj' 
(B) The following analysis of the thermodynamic and 
transport data also shows that the /? and 7 states are not 
appropriate in explaining their T-dependences. Then we 
are left with the only two states: the bipolar and pla- 
nar states in Table 1. The former is non- unitary while 
the latter is unitary. Both have no spontaneous moment 
averaged over the Fermi surface. We proceed in the fol- 
lowing sections with these states to calculate physical 
quantities. In Appendix we provide a proof that these 
particular states are local minimum and thus stable so- 
lutions of the general Ginzburg-Landau type free energy 
functional under a certain parameter space spanned by 
Pi ~ Ps because it is known that the bipolar and pla- 
nar states withpUaree-component orbital part are saddle 
point solutions .lilP Then we assume that either bipolar 
state or planar state is stabilized by fulfilling the above 
conditions (see Appendix). 

2.2 Gap function and gap topology 

The spin triplet pairing function is described by a 2 x 2 
form 



A„;3(k) = 



-dxik) +idy{k) 
4(k) 



4(k) 

da:{\s.) +idy{\i.) 



(1) 



in terms of the d(k) — (dxCk), dy(k), dz(k)) vector which 
fully characterizes each state ((ii(k) is a complex num- 
ber). The gap function is given by 



A(k) oc V|d(k)|2±|d(k)xd*(k)|. 



(2) 



If d(k) is non-unitary, the gap function has two branches, 
each corresponding to the spin up and down pairs. If 
d(k) is unitary, the gap function has only one degenerate 
branch given by 



A(k) « |d(k)|. 



(3) 



We caJl the bipolar state belonging to a non-unitary 
stateEa) when d(k) = (A:r(k), iAy(k), 0) and the pla- 
nar state belonging to a unitary state when d(k) — 
(XxCk), Aj,(k), 0) where the two kinds of the orbital func- 
tions Aa;(k), and A2;(k) are real number functions. The 
orbital function (Aj;(k), iAj,(k)) must belong to two- 
dimeufiional irreducible representations in Dgh, namely, 
Ei„,E3' or E2u- Since there must be a line node on the 
equator of the Fermi surface, we are bound to choose the 
orbital function in E2„ given by 



\x{k) ^ k^kl - kl)Gtotal{k) 

Ay(k) = kz2kxkyGtotaii'k) 
Gtotaii^) = {I + a2kl + a4kt) 



(4) 
(5) 



{i+b2{ki+kif+b,{ki+kir}. (6) 



These are supplemented with a totally symmetric func- 
tion Gtotai(k.) parameterized by 02, 04, 62, and &4. 

The actual gap topology for the bipolar state is easily 
found from Eq. (2); in the polar coordinates (0,0) on 
the Fermi sphere 



^A2(k) + A2(k)±2A.(k)A,(k) 
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I COS 6*1 sin^ Q\J\ isiii'; 



(7) 



indicating that (1) four line nodes parallel to the fc^-axis, 
(2) one line node on the equator in the absence of a2, 04, 
62, and &4, which are chosen to tune the fine details of 
the nodal topology. (3) Since in the A and C phases one 
of the two components vanishes, there are one line node 
on the equator and two line nodes parallel to the /c^-axis 
whose positions are different for the A and C piiases. 

In the unitary planar state the gap topologycHi' is de- 
termined by 



A2(k)+A2(k) = |cos6'|sin2 6', 



(8) 



in the absence of a^-, 04, 621 a-nd 64. This shows (1) the 
two quadratic point nodes at the poles (d = 0,7r), (2) a 
linear line node on the equator {9 = 7r/2), and (3) in the 
A and C phases there are either two point nodes or one 
line node. 

The T-dependence of the gap function A j. (k, T) for 
the bipolar state is given by 



where 



A2±(k,T) = A2(k,T) + A2(k,T) 
±2A,(k,T)Aj,(k,r) 



A.(k,T) = A(-)( — )A,(k) 



- a(3')| 



T 



A,(k,T) = A^«^( — )A,(k). 



(9) 

(10) 
(11) 



The order parameter amplitudes Aq ^(T) and Aq ■^(T), 
where due to the T^ splitting {Tc2/Tci = 0.9) each ampli- 
tude appears either at T = Td or T = Tc2, are evaluated 
by solving the BCS type gap equations. 

Similarly, in the planar state case the gap function is 



A2(k,r) = A2(k,r)-f A2(k,r) 



(12) 



instead of Eq. (9) in the bipolar state. The order param- 
eter amplitudes Aq -^(T) and Aq ^(T) are given by Eqs. 
(10) and (11). 

§3. Thermalconductivity and ultrasound atten- 
uation 

In order to identify the proposed two states, and to see 
which is better, we calculate their transport properties: 
the thermalconductivity and ultrasound attenuation co- 
efficient taking into account resonance scattering effects 
in the unitarity limit where the phase shift ^0 = 7r/2. 
The formulation foe. these quantities just follows those 
by Hirschfeld et a£3 and by Norman and Hirschfeldu^ 
who extend the former to considering inelastic scattering 
effects important in high temperatures near Tc. We are 
interested in the clean limit behavior here. 

The thermalconductivity Hi{T) for the heat current 
along the i-direction is given by 

^i^{T) _ 9 



i^n{Tci) 



27r2r, 

(fcfRe^ 



du^{^fsech^^K,{uj,T){l2.) 



|A(k,T)|2 



(Re 



(14) 



where k,n{Tci) is the thermalconductivity of the normal 
state at T = Td. The average over the Fermi surface 
is denoted as (•••). We assume the Fermi sphere in the 
following, ki is the unit vector for the i-direction. 

The ultrasound attenuation coefficient a^j (T) for the 
sound wave with the polarization parallel to the i- 
direction and propagating along the j-direction is given 

by 



Kr) 



aN{Tci) 



1 
1 



dujsech^—A,j{uj,T) 



(15) 



(n?.-ReV^2_|A(k,T)|2) (16) 



"^j — KiKj 






(17) 



V'^2-|A(k,T)|2' 



where a^lTd) is that in the normal state at T^i- These 
formulae are standard ones except that (1) when treating 
non-unitary state, we calculate these quantities for the 
up and down spin pairs independently, and then sum over 
the two contributions. This is permissible because these 
physical quantities heat and sound conserve the spin. 
(2) According to Norman and HirschfeldcJ the effects 
of inelastic scattering in the normal state are taken into 
account by scaling the results of Eqs. (13), and (15) by 
a factor(l -t-/)T^^i)/(l + bT'^) with b being an adjustable 
parameter.c3 



We show in Fig. 1 an example of the calculated results 
for the bipolar state (a) and the planar state-(ib) together 
with the experimental data by Lussier et al.Ell' Although 
we do not exhaust all combinations of various parameters 
(02, 04, 62, 64, b), both fittings are equally good. In Fig. 2 
the T-dependence of the thermalconductivity anisotropy 
Kc/i^b for both states arejXpmpared with the experimen- 
tal data by Lussier et alES and Suderow et aLES* It is 
rather difficult to decide which is better. 

The transverse sound attenuation coefficients for the 
two kinds of the sound propagations are plotted in Fig. 3 
(the bipolar state (a) and the planar state (b)) with the 
same set of the parameter values as in the thermalcon- 
dtictivity and compared with the data by Shivaram et 
al.c^ It is seen from Fig. 3 that both theoretical curves 
for the bipolar and planar states can explain the overall 
T-behaviors. Also notice that the theoretical curves show 
a knee structure at the second transition Tc2/Tc]|-^ 0.9. 
This anomalvjs clearly observed by Taillefer et alEZP and 
EUman et aLcS and not so clear in Shivaram et al's data 
cited here.Ea' We emphasize that although we did not 
hunt all the parameter space, both physical quantities 
K and a for two directions are consistently explained by 
the same set of the parameter values. 



§4. Phase transitions 

Combining the orbital symmetry and the previous 
spin symmetry, we have completely determined the pair- 
ing symmetry realized in UPta. We are now in po- 
sition to explain the phase diagram observed experi- 
mentally although the essentially same explanation is 
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Fig. 1. The T-dependence of the thermalconductivities Ki{T)/T 
normahzed by its value at T = Td for the fe-axis and c-axis. 
(a) The bipolar state (02 = 0.2, 04 = 3.8, 62 = 0.4, 64 = 0, 
and b = 5.0/i^^), and (b) the planar state (02 = 4.4, 04 = 0, 
62 = 1-8, 64 fhO, and 6 = ll.O/A'^). The data are taken from 
Lussier at aLrT 



given before.EH' The three basic assumptions are aheady 
adopted in this paper: (1) The effective spin-orbit cou- 
phng feh by a Cooper pair is weak. (2) The spin of 
the pairing is triple-degenerate. (3) The antiferroniag- 
netic fluctuations characterized by the triple-Q vectors 
with Qi = (i,0, 0) and its equivalent positions Q2 and 
Q3 in reciprocal space (a* = b* = -^y^, c* = ^) in the 
hexagonal planeEj' are responsible for breaking the triple 
degeneracy of the spin. 

These assumptions lead to the following GL free en- 
ergy 



F = F(2)+f(4) (18) 

^(2)= Y. «o(r-ri)Kf (19) 

j—a^b.c 

where F^'^' is given in Appendix for the bipolar state; 
Eq. (A-9) and planar state; Eq. (A-10). These par- 
ticular quartic terms ensure that the bipolar or planar 
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Fig. 2. The T-dependence of the anisotropy Kc/ki, normalized 
at T = Tci for the bipolar (dashed line) and the planar (solid 
line). The parameter values a#« the same as in Fig. 1. The dp*«, 
are taken from Lussier at alcj^ (circles) and Suderow et a1E3 
(diamonds). 



states is stable in the ground state identified as the B 
phase and that the two component spin vectors are or- 
thogonal to each other. Here we have taken the order 
parameter as db = 6Aa;(k), dc = cAj,(k) and da = aAj,(k) 
for the planar state and as db = bX^Ck), dc = icAy(k) 
and da — idXyCk) for the bipolar state. The three tran- 
sition temperatures T° > T^ > T^ which are originally 
degenerate are split into the two groups: T^ and {T^ and 
T^} by the symmetry breaking field as mentioned. The 
latter two are further assumed to be slightly different: 
T° < T^ due to the small uniaxial anisotropy of the sys- 
tem. At T = T^!, which is identified as the upper critical 
temperature Td = 0.58K the A phase characterized Jav 
da = 0,db T^ and dc = appears first (see Fig. 4).cj) 
Then at a lower temperature, T = T^ which is identified 
as Tc2 = O.SSiiT the second order transition from the A 
phase to the B phase which is characterized by either the 
non-unitary bipolar state or the unitary planar state. It 
can be proved within the above free energy that third 
transition at T — T^, which is designed to situate a few 
mK below T = T^ never realized at zero field. 

In order to discuss the phase diagrams for external 
field H applied parallel to and perpendicular to the_c- 
axis. We must take into account the vortex structures' 
based on the above free energy functional by adding the 
terms describing the spatially varied order parameter: 



-^ ^ grad 1" ^ hulk 



Fgrad^ J2 {Ki{\D,d,\^ + \Dyd,\') 
j—a,b,c 

+Ki\D,J,\'} 
Fbuik= Y. MT-W\dj\^ 



(20) 



(21) 



j—a,b,c 



^AxplH-dp + i^W. (22) 



The gradient term Fgrad describes the spatial variation of 
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Fig. 3. The T-dependence of the ultrasound attenuation coeffi- 
cients axy{T) and a^^iT) normalized at T = T^i- (a) The 
bipolar state, and (b) the planar state. The parameter values 
are-the same as in Fig. 1. The data are taken from Shivaram at 
al.Ei) 



the order parameter under external field with the gauge 
invariant derivative Dj = —ihdj — —Aj (A is the vec- 
tor potential). Here there exist two kinds {Kf and K2 
{j = a,b and c)) of the gradient term associated with 
hexagonal symmetry where K^ {K2) can differ from Kf 
and Kf {K2 and K2) because of the symmetry break- 
ing AF fluctuation, that is, K^ ^ Ki - CM^,;^^ = 
K2 - C'M2, K'^ ^ Ki+ CM^.K!^ = K2 + C'M2, Kl = 
K[+CyP,Kii = K'2+C,'M.^. Here M^ is the amplitude of 
the AF fluctuation polarized along the 6-axis. It should 
be noted that there is no so-called gradient mixing term, 
which washes out the desired tetra-critical point, in the 
present spin scenario. jTiiia is quite different from that 
in the orbital scenariosQ'B'o' where the gradient mixing 
is inevitable. 



The resulting phases in the H ys T shown in Fig. 4 
are fully characterized by specifying the d vector in each 
phase, that is, its direction, the number of the compo- 



H II a,b 




Tc£ Tci 




Fig. 4. Schematic phase diagrams in _ff vs T for H || a,b and 
H \\ c where the direction and the components of the d vector 
for each phase are indicated. Hgc is the second phase transition 
between the B and C phases. Hrot is the rotation field of the d 
vector 
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Fig. 5. Schematic figure of the spin rotation under a weak H 
field. 



nents and the associated orbital parts. The three vectors 
in the spin-triplet state are now denoted by the real vec- 
tors a, b and c referring to the Cartesian coordinate in 
the hexagonal crystal. The high T and low (high) field 
A(C) phase is described by a single component d vector, 
while the B phase in the low T and low H is a non-unitary 
bipolar state or unitary planar state characterized by 
two-component d vector. Because of the Axp term in 
the above free energy Eq.(22), which tends to align the 
d vector perpendicular to H, the B phase for H || c is 
further subdivided into bAa;(k) -I- (i)cAy (k) in low H and 
bAa;(k) -I- (i)aAj,(k) in high H. As H{\\ c) increases the d 
vector rotates, corresponding to the Knight shift behav- 
ioiO at H ~2kG. The d vector rotation is schematically 
shown in Fig. 5 where one of the d vector components 
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rotates upon increasing _ff(|| c) while keeping the other 
component fixed to the 6-axis. 

§5. Summary and conclusions 

Before going to a summary and conclusions, we give 
additional discussions to clear our standpoints ar 
answer possible objections raised in several places! 

5.1 Supplemental discussions 

The double transition : The origin of the Tc-splitting 
comes from the crystalline hexagonal symmetry break- 
ing due to either incommensurte lattice modulations, or 
stacking faults along the c-axis occurring from a high 
temperature, which are jijell correlated with the exis- 
tence of the Tc-splittingE2l^ This leads to the different 
spectral weights of the otherwise triply degenerate anti- 
ferromagnetic fluctuations whose polarization points to 
the a* (II 5)-axis, thus the originally three equivalent a*- 
axes become inequivalent from a high T . This inequiv- 
alence is mdecd observed in Ton et al's Knight shift ex- 
periment. £3 The spectral weight difference of the anti- 
ferromagnetic fluctuations breaks the degeneracy of the 
multinlfi rnrnrinripnt; order parameters. The growing ev- 
idcnccEj'OoEj-' shows that the antiferromagnetism is 
not true static and long-range order, which is unable to 
work as a symmetry-breaker as postulated by the Eig 
and E2m scenarios. This is gumain reason why we aban- 
doned our original scenaricEa) Moreover, the expected 
antiferromagnetic moment rotation upon applying the 
basal plane magnetic field is not observed by the neutron 
experiment. EJ' This makes these scenarios more diffi- 
cult because the basal plane phase diagram is anisotropic 
which is not seen experimentally. As for the main pairing 
mechanism we think the observed ferromagnetic fluctua- 
tiona23^ for stabilizing the triplet channel as in superfluid 
■^He while the antiferromagnetic fluctuations mentioned 
above are merely to split the transition temperature with 
otherwise degenerate state. Indeed the latter fluctuation 
intensity is known tojificrease when the system under- 
goes superconductingJSj' revealing it workable as a sym- 
metry breaker. 

The spin rotation: It was remarkable to observe that 
the spin susceptibility begins to decrease as H(\\ c) in- 
creases, implying that the d- vector rotates at around 2kG 
within the B phase upon responding to external field to 
lower the Zeeman energy. This immediately leads us that 
(l)the effective spin-orbit coupling felt by a Gooper pair 
is weak contrary to what has been claimed.cHP (2)The 
spin component of the order parameter must be nn. 
tiple, negating not only the spin-singlet scenario Ei^j 
but also the single spin component scenario such as 
which also assumes the strong spin-orbit coupling.! 
Our scenario is able to explain this remarkable observa- 
tion as explained before. 

The Hc2 crossover phenomenon: The upper critical 
field H^ for the basal plane exceeds that for the c-axis at 
low T. The T-dependence of H^ exhibits no saturation 
towards T — > 0. This crosjsaver behavior of H^ and iJJ,2 
is often taken as evidenceliil' of the absence of the Pauli 
limiting of H^. We think, instead, that the anomalous 



H^ curve indicates a new phase above the pjC phase. The 
recent Knight shift experiment by Tou et aSi3 shows that 
X{H) begins to decrease as H{\\ a) increases at around 
lOkG in the middle of the C phase. 

5. 2 Conclusions 

We have demonstrated that the proposed orbital 
part combined by the previously determined spin part, 
namely the non-unitary bipolar state and unitary pla- 
nar state gives rise to a reasonable explanation for the 
overall T-dependences of the thermalconductivity and 
transverse ultrasound attenuation both for the c-axis and 
the basal plane. These states can also explain the ma- 
jor previously known thermodynamic and transport ex- 
periments, including the specific heat (^ T"^), the nu- 
clear relaxation rate (~ T^) and the longitudinal ultra- 
sound attenuation for both directions (~ T) , all of which 
merely point to the existence of a line node since the 
former two are directionally averaged quantities. The 
directioi^resolved anisotropic penetration depth probed 
by ^SRcS' for two direction (|| c and _L c) is interest- 
ing to interpret in terms of these states, which belongs 
to a future study. Other measurements related to sur- 
face sensitive_technique such as point conlact tunneling 
spectroscopyo* and Joscphson tunnelingpj^ which yield 
valuable information of each phase A, B and C phases 
remain to be done although we believe that our bipolar 
or planar states are basically able to understand these 
data. 

It is rather difficult to distinguish between the bipo- 
lar and planar states both theoretically and experimen- 
tally. The existing Knight shift experiments are consis- 
tent completely with both states. There is no distinction 
between them. As for the orbital part related to the gap 
topology in the B phase the bipolar state has four line 
nodes parallel to the c-axis while the planar state has 
two quadratic point nodes at the poles in addition to 
a line node on the equator in common. Thus we expect 
that there exists the in-plane four-fold anisotropy for the 
bipolar state while there is no in-plane anisotropy for the 
planar state. We warn, however, that this anisotropy 
may not be so strong to be able to discriminate them 
because of the background hexagonal symmetry. But 
this kind of experiments, including the transverse sound 
attenuation for various in-plane polarizations under the 
sound propagations kept along the c-axis are worth try- 
ing. In connection with the iJc2 crossover phenomenon 
mentioned above the new phase transition which ap- 
pears when H{\\ a) ~ lOkG should be detected by other 
phase transition sensitive experiments such as elastic 
anomalycS* It is also, interesting to further investigate 
the vortex structurecS in terms of the identified state. 

Appendix: Stability conditions of the bipolar 
and planar states 

We start with the general free energy functional valid 
up to the quartic termstZP invariant under spatial and 
spin rotations and a guage transformation given by 
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/(4)=/3itr|Ai|2 + /32(tr|AAt|)2 
+/33tr[(AtA)(AtA)*] 

+PMAA^? + PMiAA^)iAA^)*] (A-3) 

where A^i is a complex 3x3 matrix. /i(= 1,2,3) and 
i{— x,y,z) denote the spin and orbital indices respec- 
tively. A and A^ are the transpose and Hermitian conju- 
gate of ^. In the case of the p-wave pairing such as super- 
fluid He"^, A relates to the d- vector as d(k) = A-k. When 
the orbital component is confined to just two, whose ba- 
sis functions arc XxO<-) and Aj,(k), we can reduce the 
d- vector as d(k) — A^^xO^) + ^y^yO^)- 

By writing the d-vector as A^ = d+if3, and Ay = ^+id 
where a, /3, 7 and S are real vectors, we find from Eq. 
(A-3) 

f^'^'> =/3i{(a-a-/3-/3 + 7-7-(5-(5)2 

+4{d-p + ^-S)} 
+l32id ■ a + (3 ■ (3 + J ■ J + S ■ Sf 
+/33{{^ -d + p-pf + ii-j + S-Sf 

4-2(a-7 + /3-(5)2-2(/3-7-a-5)2} 
+f34{{d -d + p-pf + ij-j + S-Sf 

-f-2(a • 7 + /3 • (5)2 4- 2(/3 • 7 - 5 • 5)2} 
+f35{{d • a - /3 • /3)2 + (7 . 7 - 5 • 5)2 

+4{d ■ pf + 4(7 • 5)2 + 2(a • 7 - /3 • 5)2 

+20-^ + d-5f}. (A-4) 

In order to know the stability of the bipolar and planar 
states, let us consider small deviations from the respec- 
tive states: A^i = A°j-|-a(y^, i)+ib{^, i) with i = x,y. For 
the bipolar state /bipolar is written up to the quadratic 
order in small deviations as 

SfMpoiar = /3i{8a2(l, x) + 2(a(2, y) + b{l,x))^} 
-|-/33{4a2(l,a;) + (a(l,y) + 6(2,a;))2 

-(a(2,.T) + 6(l,2;))2} 
+/34{4a2(l, x) + (a(l, y) + 5(2, x)f 

+ {a{2,x) + bil,y)f} 
+/35{4a2(l,a;)-(a(l,y) + 5(2,a;))2 

+ {a{2,x) + bil,y)f 

-2a2(3,y)- 262(3,2;)}. (A-5) 

This leads to the stability conditions for the bipolar 
state: 

/3i > 0, /?5 < 0, 2/3i + /33 + /34 + /35 > 

P3 + I34- P5 > 0, -P3 + /34 + /35 > 0. (A.6) 

For the planar state Sfpianar is calculated as 

Sfplanar = /3i{-2(6(l, x) - 6(2, y)f - 462(1, y) 

-462(2, .t)- 462(3, a;)- 462(3, y)} 



triplet pairing UPtg 7 

+/33{4a2(l,.T) + (a(l,y) + 6(2,a;))2 

-(6(l,y)-6(2,x))2} 
+p4{^a^{l,x) + (a(l, y) + a(2, x)f 

+ (6(l,y)-6(2,x))2} 
+A5{4a2(l,a;) + (a(l,y) + a(2,a;))2 

-(6(l,y)-6(2,x))2 

-262(3,x)- 262(3, y)}. (A-7) 

The stability conditions for the planar state is given by 

/3i < 0, /35 < 0, 

/33 + /34 + /35 > 

-/33 + /34 - /35 > 0, (A-8) 

The resulting quartic terms for each state is rewritten 



as 



/Soia. = (/3l + 2/5345)(5-^-'?-^T 

+ {132 + ^P345){d ■ d + S ■ Sf 
+2(-/33 + /345)(a-^)' 

4La. = (/3i2 + ^/5345)(a-5 + 7-7)' 

1 
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+ ^/3345{(a • a - 7 • 7)2 + 4(5 • 7)'}A-10) 

where Pijk--- = A + /3j + /3fe + • • •• The minimization is 
found to require from Eqs. (A-9) and (A- 10) that a ± 5 
and \d\ — \5\ for the bipolar state and a _L 7 and \d\ — \^\ 
for the planar state. Namely, the spin component vectors 
for each orbital component are orthogonal to each other. 
The general structures of the 3x3 matrix A under a fixed 
a vector parallel to the 1-direction arc now written as 
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bipolar 



ai 
i52 
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(A-11) 



ai 

Ar,lar,.ar = I 72 ) (A- 12) 

73 

where ai, 52, 53, 72 and 73 are real numbers. In the 
main text the spin indices /i~l, 2, 3 correspond to the 6, 
c and a-axes respectively. 
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